Quantum Hall effects of graphene with multi orbitals: 
Topological numbers, Boltzmann conductance and Semi-classical quantization 
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Hall conductance a xy as the Chern numbers of the Berry connection in the magnetic Brillouin 
zone is calculated for a realistic multi band tight-band model of graphene with non-orthogonal basis. 
It is confirmed that the envelope of a xy coincides with a semi-classical result when magnetic field is 
sufficiently small. The Hall resistivity p xy from the weak-field Boltzmann theory also explains the 
overall behaviour of the a xy if the Fermi surface is composed of a single energy band. The plateaux 
of a xy are explained from semi- classical quantization and necessary modification is proposed for the 
Dirac fermion regimes. 
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The quantum Hall effect (QHE) has been one of the 
important subjects in condensed matter physics for sev- 
eral decades. Wide range of experimental and theoretical 
studies have been devoted to this subject. In the early 
stage of these efforts, the topological aspects ofquantiza- 
tion have been firmly established^, H i, S, M, 0- The 
topological formulation of quantized Hall conductance Q 
has been utilized numerically to explore exotic nature of 
electrons on two-dimensional lattice models. 

Recently, electrons on a honeycomb lattice has at- 
tracted theoretical attention because energy bands on 
this lattice are described by massless Dirac disper- 
sion near the center of bands. The dispersion results 
in anomalous QHE a xy — (2N + l)e 2 /h (N: inte- 
ger) per spinj^, [l(|. Such anomalous QHE was ob- 
served experimentally in graphene, which is monolayer 
graphite [111, [13]. The topological aspect of QHE in 
graphene has been studied by using a single-orbital tight- 
binding model on a honeycomb lattice pi 



b| . It has 



been shown that anomalous QHE persists up to the van 
Hove singularities [3|. At the singularities, a xy jumps 
discontinuously and normal QHE a xy = Ne 2 /h is recov- 
ered around band edge regions. These results suggest 
that bands structure affects QHE both qualitatively and 
quantitatively. Thus, it is an interesting question to in- 
vestigate the relation between band structure and a xy . 

In this letter, we study relations between band struc- 
ture and quantized a xy by using a realistic tight-binding 
model of graphene as an example. Up to now, numerical 
calculations of quantized a xy have been rather limited to 
simple tight-binding models and realistic band structure 
was beyond scopes of these studies. Another purpose 
of this study is to demonstrate that the topological for- 
mulation of a xy [3, @| can be applicable to realistic band 
structures. 

We use a tight-binding model for graphene with s and 
p orbitals on carbon atoms [Hj]. The spin degeneracy is 



ignored for simplicity. Each orbital is labeled by a 
site index i and an orbital index £. Transfer and 
overlap integrals Si£,jf between nearest-neighbor atoms 
are considered by a Slater-Koster approximation. The on- 
site energy e^, tit^ fi a nd Sif.if are obtained from the 
parameters in Refs. [l5| and LLfj. With eight orbitals in a 
unit cell, energy bands dispersion is calculated as shown 
in Fig. 1. In actual graphene, electrons are occupied 
up to e = 0. Around this energy, band dispersions are 
well approximated by the massless Dirac cones, which 
is the origin of several interesting physical properties 
including anomalous quantization of Hall conductance. 
The massless Dirac dispersions are also realized around 
s ~ — 14eV. In addition, several van Hove singularities 
are clearly seen in the total density of states (DOS). 
The uniform magnetic field is introduced as 
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so that they satisfy 

Ecloscd loop % = 27r(flux quanta in the loop). 

Under the uniform magnetic field, original transla- 
tional symmetry is broken. The system is then char- 
acterized by the magnetic flux <f> = BCl/ipo where SI is 
the area of a unit cell and (pa — hc/e is the quantized 
magnetic flux. When is a rational number p/q, we can 
define magnetic unit cell whose area is q times larger than 
that of the original cell. Then, extended Bloch theorem 
holds for the magnetic unit cell and eigenstates are la- 
beled by a wave- number k. Thus, an energy band splits 
into q sub-bands by the uniform field. If the chemical 
potential is located within an energy gap, quantized Hall 
conductance a xy is given as 



-c F ([x), 



(2) 



where cf is a topological integer called Chern number de- 
fined for the filled bands. It is calculated by discretizing 
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FIG. 1: (Color online) The lower panel shows the calculated 
Chern number {cr xy ) for <f> = 1/53 as a function of the chem- 
ical potential. The upper and middle panels present energy 
bands and total density of states when magnetic field is ab- 
sent. Vertical dashed lines indicate positions of van Hove 
singularities. 



the Brillouin zone into mesh {ki} as [HID, Q 

Mm) = ^E F ( fc *) 

t 

F(k) = krgu x {k)u v {k + Ak x )[u x (k + Ak y )u y (k)]* 
u„(fc) = det[*(fc)] t *(fc + Afe M ) 
*(fc) = (IVi(fc)),-- - ,|lM*)» 

where Afc„ is a discretized momentum along v direction, 
\ip m (k)) is a Bloch state of the band index m and m = M 
is the highest occupied energy band below the chemical 
potential fi. (See Ref.[8( for the detail). 

This formulation was successfully applied to quantum 
Hall effect on a disorder system [171 and graphenel4|. 
When the field <f> = p/q is sufficiently small, p sub-bands 
are grouped together in general. These grouped bands 
correspond to a Landau band of nearly free electrons. 
Hereafter, we treat the field <fi — 1/q, which enables us 
to consider the weak-field limit easily. 

In Fig. [TJ we show the calculated Chern number cf 
with tf> = 1/53 as a function of the chemical potential. 
The Cf of adjacent gaps are connected by a straight line 
as an eye guide. In some energy regions, several Lan- 
dau bands are grouped together with very small energy 
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FIG. 2: (Color online) Comparison of the Hall conductance 
obtained from Chern number with semi-classical formulations. 
The solid line shows the a xy obtained from Chern number 
and open circles indicate that from semi-classical theory with 



clean limit. The dashed line is 
limit of Boltzmann theory. 



-p xy obtained from weak-field 



gaps. In such case, we found that Chern numbers corre- 
sponding to these small gaps are sensitive to the choice 
of (j). Thus, they are not physically relevant and we omit 
them from Fig. [1] The plateaux are most clearly seen 
at Dirac fermion regions because of large energy gaps 
between Landau bands. The cf increases 1, 2, or 3 be- 
tween adjacent (well-developed) energy gaps. If these 
steps are smeared, coarse grained a xy behaves as a con- 
tinuous function except for several energies where the cf 
changes with large numbers. 

We found that the large discontinuities in Cf can be 
attributed to topological changes of Fermi surface when 
the magnetic field is absent. For example, the discon- 
tinuity at e » — 15.7eV is caused by a transformation 
of Fermi surface from an electron-pocket around Y point 
to two hole-pockets around K points (see Fig. [3] (a) and 
(b))- 

Such topological changes are also recognized as a van 
Hove singularity in DOS. Other discontinuities also cor- 
respond to similar topological changes of Fermi surface. 

We compare global behaviour of a xy with a semi- 
classical theory. When Fermi surface consists of closed 
curves, the Hall conductance a xy for clean-limit can be 



calculated as 18, 1 
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within semi-classical approximation. Here, Si (p) is the 
directed area enclosed by i-th segment of Fermi surface. 
It is defined as a positive (negative) value when the area 
enclosed by the segment is filled (empty). Above equa- 
tion can be casted to the form: 
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where J7bz = (27r) 2 /f7 is the area of Brillouin zone. We 
found that c sc completely describes the envelope behav- 
ior of topological number cf for weak field limit as shown 
in Fig. [21 While this result may be naively expected as 
we compare same physical quantity a xy , it is highly non- 
trivial from the fully quantum mechanical formulation 
of Hall conductance. We have succeeded in demonstrat- 
ing numerically that quantized Hall conductance agrees 
with semi-classical limit when the applied field is suffi- 
ciently small. Conversely, this agreement indicates that 
the formulation based on Chern number is quantitatively 
applicable to calculate a xy in clean-limit for realistic en- 
ergy bands. Thus, this formula may be useful to explore 
quantum Hall effect for materials with exotic energy band 
structures. 

The topological changes of Fermi surface occur as 
transitions from electron pockets to hole pockets. At 
these points, J2i &i discontinuously decreases with £Ibz- 
Therefore, discontinuity of Chern number is </> -1 = q and 
Hall conductance jumps by (— e 2 //i)0 _1 universally. 

Another well-known semi-classical formulation of Hall 
effect is the Boltzmann theory within relaxation time 
approximation [2Ct l2ll |22|| . If the relaxation time r is ap- 
proximated by a constant, Hall resistivity p^ y for weak 
magnetic field limit does not depend on r and is obtained 
asl'211 
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where 

riu = (vf), r] xyz = -(vB^Vg^-Vy^Vy}, v = 
Vfe£fe is a Fermi velocity and (• • • ) implies an average over 
the Fermi surface. For free electron dispersion £k = 
the Boltzmann theory gives p® y = -^B, which coincides 
with clean-limit o & £ y = —nec/B as p^ y = —l/a^ y . There- 
fore, it would be interesting to verify whether such cor- 
respondence holds for general band structures. In Fig. [21 
we compare p^ y with Chern numbers. We found that 
Boltzmann theory generally agrees with Chern numbers 
when the Fermi surface is essentially generated by a sin- 
gle energy band and the chemical potential is not close to 
van Hove singularities. When multiple bands contribute 
to Fermi surface, interband effects become important and 
two limits of semi-classical theory do not coincide. In 
such case, the a xy would sensitively depend on the de- 
tails of scatterings or relaxations. 

Next, we examine the quantized plateaux of a xy in de- 
tail. Fig. [3] shows enlarged figure of cp for <f> = 1/200 at 
several energy regions. At the bottom of energy bands 
e ~ — 17.8eV, where Fermi surface consists of an elec- 
tron pocket around T point, cp increases with a step of 



3eV % e % 3eV are shown in Fig. [3] (b) 



1. Therefore, the quantum Hall effect in this region re- 
sembles that of free two dimensional electron gas. 

Massless Dirac fermion regimes at 
— I3cV and 

and (c). In these regions, Fermi surface is approximated 
by two massless Dirac cones centered at K points. Here, 
two Landau bands are grouped together and the energy 
gap between them are extremely small. Well-developed 
energy gaps appear between these grouped Landau bands 
and the Chern number cf increases with 2 between ad- 
jacent grouped bands. This behaviour may be explained 
from the existence of two equivalent Fermi surface seg- 
ments. However, naive interpretation would give cp = 2n 
which differs with actual steps cf = 2n + 1. These 
plateaux have been theoretically obtained by solving the 
Dirac equation under magnetic fields flo| or using the sim- 
plified tight-binding model for carbon 7r-orbital[Ii HI]. 
In the present calculation, we included realistic band 
structure which deviates from the linear dispersion Ek ~ 
±?j|fc|. The Chern numbers do not change by such de- 
tails because topological number is robust as far as en- 
ergy gaps do not close. However, the positions of steps 
are influenced by the band structure. It is most clearly 
seen in Fig. [3] (c) as a violation of the electron-hole sym- 
metry around e — 0. Namely, the widths of plateaux for 
e < are narrower than those for e > 0. The quantitative 
explanation of such behaviour will be given below. 

Near e ~ 18eV, Fermi surface consists of three electron 
pockets around M points. We found that three Landau 
bands are grouped together and Chern number can be 
written as cp = 3n. More interesting behavior is real- 
ized in other energy regions. For example, Fermi surface 
around e ~ 20eV is composed of a hole pocket around T 
point and two hole pockets around K point. In this case, 
cf shows two differnet steps with increments of I and 2 as 
presented in Fig. [3] (e) . This behavior can be understood 
by an assumption that cp increases I (2) when the semi- 
classical quantization condition to Fermi surface around 
T (K) points is satisfied. 

Let us perform quantitative comparison with the semi- 
classical quantization. Onsager's semi-classical quantiza- 
tion condition 23j can be written as 
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We assume that Fermi surface segments individually 
quantize by above condition. When the chemical poten- 
tial is located between quantized levels 



Si (u) 

n - I + 7 < — — < n + 7, 
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we further assume that Hall conductance a xy is given 

2 

by a quantized value — \n. Summing up a xy from all 
Fermi surface segments, we arrive an expression for Hall 
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FIG. 3: (Color online) Chern numbers for <f> = 1/200 at several energy regions. The upper panels show the corresponding 
Fermi surface calculated at the energies indicated by vertical lines. The shaded regions are filled by electrons. The dashed 
lines indicate the semi-classical estimation c^ q denned by Eq. [5] For (b) and (c), the estimations by Eq. [5] are also plotted but 
invisible because they trace almost same points with solid lines. 



conductance 
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where [_" " J denotes the largest integer not exceeding 
its argument. If we choose 7=5, c scq explains quantized 
Hall conductance in several energy regions. In Fig. [31 the 
c scq is plotted as dashed lines. They completely agree 
with cf for Fig. [3] (a) and (d). When multiple types of 
Fermi surface segments exist, the semi-classical quanti- 
zation deviates from calculated cf- It may be caused by 
multi-band effects. 

At massless Dirac fermion regions, the c scq fails to ex- 
plain cf because semi-classical quantization near mass- 
less Dirac cones is different [241] . The existence of topo- 
logical phases influences the quantization 25, 2f| 13] • Re- 
cently, P. Gosselin et a?[28[ have shown that semiclassi- 
cal quantization condition of massive Dirac fermion is 
described by 7 = and 1. Extrapolating this result to 
massless limit, we obtain a modified expression of c scq for 
Dirac fermion: 
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where So(fJ') is the area enclosed by a Fermi surface 
segment around K point. This equation is also plotted 
in Fig. [3j (b) and (c) but almost invisible because it com- 
pletely traces the Chern number cp. Thus, we found 
an expression of the Hall conductance for Dirac fermion 
regime. It completely predicts the positions of quantized 



plateaux and explains asymmetry around e = 0. If bias 
voltage toward graphene can be sensitively controlled, 
the prediction by Eq. [9] may be compared with experi- 
mentally observed plateau structure. 

In summary, we demonstrated that quantized Hall con- 
ductance of general band structure can be practically 
calculated from topological numbers defined to energy 
bands under magnetic fields. Using graphene as an exam- 
ple of multi-band models, it was shown that envelope of 
quantized Hall conductance coincides with semi-classical 
expressions when magnetic field is sufficiently weak. The 
semi-classical quantization can predict the positions of 
plateaux for simple bands. For Dirac fermion regions, 
modified expression was found to predict the plateaux 
completely. This expression may be utilized for analysis 
of experimental data. 
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